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Abstract 



This study investigates the potential of the TESLA Linear Collider for measuring reso- 
nance parameters of Higgs bosons beyond the mass range studied so far. 

The analysis is based on the reconstruction of events from the Higgsstrahlung process 
e"'"e~ HZ. It is shown that the total width Fh, the mass mpj and the event rate can be 
measured from the mass spectrum in a model independent fit. Also, the branching ratios 
BRh-^ww and BRh^zz can be measured, assuming these are the only relevant Higgs decay 
modes. 

The simulation includes realistic detector effects and all relevant Standard Model back- 
ground processes. Results are given for mn = 200 — 320 GeV assuming JC = 500 fb^^ 
integrated luminosity at collision energies of ^/s = 500 GeV. 



1 Introduction 

During the past years many simulations have been performed to investigate the prospects of 
measuring Higgs boson properties at future Linear Colliders ^ [21 Ol- The main focus was 
set to Higgs masses mn below the WW-threshold which is the mass region prefered by recent 
electroweak data. 

In high energy e''~e~ collisions Higgs bosons can be produced in two dominant production pro- 
cesses: Higgsstrahlung e"'"e~ — > HZ and WW-fusion e"'"e~ Hi/eZ^e, see Fig.^ For lower collision 
energies y/s, Higgsstrahlung with cjhz ~ 1/s dominates, while WW-fusion becomes the major 
production mode due its (JuucUc ~ logs cross-section rise if -y/i is large compared to mn- The 
cross sections of both processes for the considered range of Higgs masses are plotted in Fig. [21 
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Figure 1: Higgs production diagrams: Higgsstrahlung (left) and WW -fusion (right). 



Higgs bosons couple to mass and therefore decay in general to the heaviest particles possible. In 
the Standard Model (SM) and most of its extensions, the total Higgs decay width is expected 
to be very small for Higgs masses below the WW-threshold. Fig. [31 shows the SM prediction. 
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Only if the Higgs width is as large as few GeV can it be determined from the observed Higgs 
lineshape. This is not possible for smaller widths due to limited detector resolution, but the 
width can only be determined indirectly via the Higgs couplings [HI E] • 
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Figure 2: Higgs production cross sections as predicted by the Standard Model. Values are 
calculated using HZHA ^4^. 

From the lineshape not only the width, but also the mass and event rate can be determined in 
a model-independent way. The signal and background processes studied are specified in Sec. [21 
Event selection is described in Sec. |31 followed by the methods of estimating detector resolution 
in Sec. IH For this, a separation of H ^ WW and H ZZ decays is necessary, which can be 
interpreted as branching fraction measurements assuming that there are no other major decay 
modes. The note continues with details on the reconstruction of Higgs resonance parameters in 
Sec. |SJ Results are summarized and discussed in Sec. IHl 




Figure 3: Standard Model prediction for the total Higgs width Ftot- Also shown are the partial 
widths r^jj^, Fww '^''^d Tzz of the dominant decay modes. Values are calculated using HDECAY 
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2 Signal and Background 



Higgs bosons with Standard Model couplings are studied. The parameter space of interest 
ranges from mji = 200 GeV to 320 GeV, where both cross section and width are large enough 
for precision measurements. 

SM Higgs bosons in the mass range under study, decay almost exclusively to pairs of massive 
gauge bosons. The final state of interest therefore is e~''e~ H (Z) — > WW/ZZ(Z). All 
successive decay modes are listed in Tab. ^ 
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WW 
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SM branching 
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1. 


2v + X 


X 
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35.08% 
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2. 


qq qq qq 


X 


X 


34.3 % 


6000 


3. 


qq qq lu 


X 




19.18% 


3355 


4. 


ii qq qq 


X 


X 


7.84 % 


1370 


5. 


iiqqiu 


X 




2.94 % 


515 


6. 


Uii qq 




X 


0.63 % 


110 


7. 


uuu 




X 


0.03 % 


5 



Table 1: Signal final states and their occurrence. Numbers given are estimates using BRh-^ww = 
0.7, BRh-+zz = 0.3, BRw^^i/ = BRz^^£ = 0.3, BRw-^qq = BRz^qq = 0.7 and cjhz = 35fb as 
predicted by the SM for mn = 240 GeV. Here, £ is any charged lepton flavor and q any quark 
flavor kinematically allowed. 

Approximately one third of the HZ-events contain more than one neutrino and thus a precise 
mass reconstruction is difficult. For the rest, the fully hadronic final state is by far dominant. 
But also, this is the channel where the largest background contributions are expected (e.g. 
e+e~ — >■ tt). The same is true for the next-to-dominant decay mode (hadronic Z-decay plus 
semi-leptonic decay of the W-pair). On the contrary, the gold-plated channel with six leptons 
in the final state is too rare for precision measurements. Channel 4 with one leptonic Z-decay 
and hadronic W- or Z-pair decays is a good compromise between signal rate and background 
contamination. 

Since r-decays always include neutrinos, only i = e, n are considered to achieve the best mass 
reconstruction. So, from here on lepton only means electron or muon. 

All background processes with two charged leptons plus jets are considered. They can be clas- 
sified as follows: 

6f Six fermion proceses, e^e~ — > 2i 4q, yield events with the same final state as signal events. 
As will be shown later, this is the dominant class of background processes. 

4f Four fermion processes including ZZ-pair production, e^e~ 2i 2q. First of all, this class 
of processes is problematic because of huge cross sections. However, event topology differs 
from signal events. 

tt Top quark pair production, e'''e~ ^ tt — > bbW'''W~. Here, high energetic leptons might 
not only occur in W- but also in b-decays. Therefore, all decays W — > qq and W — > 
are considered. Again it is not event topology but large cross section which makes this 
backgroud possibly dangerous. 

Other processes (e.g. e^e^ WW qqiu) are expected to be negligible due to missing isolated 
leptons, large missing energy or low mass of the hadronic system. 
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Signal as well as background events are generated using WHiZard 1.22 [H], except tt-events for 
which PYTHIA 6.2 ,9^ is used. Both initial state radiation ISR and beamstrahlung jJUl are taken 
into account. For this analysis, no significant signal over background enhancement is expected 
for polarized beams, so the possibility of beam polarization is not studied. Cuts on fermion-pair 
invariant masses mu/mqq > lOGeV for any lepton-/quark-pair are applied on MC level to save 
CPU power. These events anyhow are far from the parameter space of interest. 

Detector response is simulated using the fast Monte Carlo SIMDET 4 JTP which parameterizes 
detector performance as described in the TESLA TDR 

3 Event Selection 

In all events, energy flow objects, which are classified as lepton (electron or muon) are searched 
for. The most energetic of these objects is combined with any other identified lepton. The pair 
with invariant mass closest to mz is selected as Z — > candidate and removed from the event. 
All other energy flow objects are forced to four jets by the Durham recombination scheme |12| . 

The distributions of the most important variables used for background suppression are displayed 
in Fig. 01 Each event is required to pass the following cuts: 

1. The event must contain at least two energy flow objects classified as electron or muon by 
the detector reconstruction, > 2. 

2. Both leptons from the Z-decay must satisfy jcos^^j < 0.99, where cosOi is the lepton's 
polar angle. 

3. All jets must satisfy |cos6'jet| < 0.95. 

4. The two lepton invariant mass must be close to the Z-mass, {mu — mz| < 5GeV. This 
reduces background events with non-resonant lepton pairs. 

5. The sum of the energy of both leptons must be significantly less than the beam energy, 
Eu < 225 GeV = 0.45 y/s. In events from Z-pair production, each Z-boson carries half of 
the event energy, so events of this type are reduced. 

6. The hadronic system must be four-jet like, 7/34 > 10~^. Here, 2/34 is the separation param- 
eter between three (small 1/34) and four (large 7/34) jets of the jet finder. This cut reduces 
most of the remaining four-fermion events with non-resonant leptons still left after the rriu 
cut. 

In addition, a kinematic fit with four constraints is applied. The aim is to improve the mass 
resolution and to group the four jets into two pairs. The constraints are as follows: 

1.-2. Conservation of transverse momentum: 

3. Conservation of energy and longitudinal momentum allowing for one initial hard photon 
in ibz-direction: 

4. Four jets form two pairs of equal mass: 
mi,2 = m3,4. 
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Figure 4: Distributions of the di-lepton mass m^i, di-lepton energy Eu, jet separation parameter 
7/34 and the ~)(^ -prohahility of the kinematic fit PROB{x^). Cuts are indicated by arrows. 
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ignal mn [GeV] 




Background 




Variable (range) 


200 


240 


280 


320 


6f 


4f 


tt 


Events / 500 fb"' 


1120 


880 


630 


410 


4340 


392 000 


240 000 


1.-3. 


A'"^, cos 


01, COs6'jet 


904 


706 


498 


322 


1803 


38 000 


90 000 


4. 


Ell 


(< 225 GeV) 


897 


705 


497 


321 


1512 


5764 


90 000 


5. 


mi 


(mz±10GeV) 


770 


600 


425 


275 


369 


1506 


750 


6. 




(> 10-3) 


745 


581 


413 


268 


343 


5 


357 


7. 


PR0B(x2) 


(> 10-6) 


585 


463 


333 


216 


271 





4 


Efficiency 


52% 


53% 


53% 


53% 





Table 2: Evolution of event rates through the selection. Main remaining background source are 
events from 6f processes — > 2£4q. Event rates in the first line include MC level cuts as 
described in the text. 
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This fit is performed for all three possible assignements of the four jets to two pairs in the fourth 
constraint. The combination yielding the best fit is used. Events with a probability of this 
below PROB{x^) < 10^^ are rejected. 

Tab. 121 shows the overall performance of the event selection. Background suppression is possible 
to S/B 1 or better depending on the Higgs mass. Events from four fermion processes and 
tt-pair production can be rejected almost completely. Signal efficiency is stable as a function of 
mn and lays above esignai > 50%. 

In the following, only three objects are considered per event: A reconstructed Z-boson (the two 
leptons) and two further objects assumed to be either a pair of W- or Z-bosons (the jet pairs). 



4 Separation of H ^ WW and H ^ ZZ 

Since it is a priori unknown which two of the three bosons originate from the Higgs decay, a 
distribution is formed which contains all three possible di-boson masses per event. It is expected 
that the correct pairing will exhibit a mass peak while the two wrong pairings will form a flat 
combinatorical background. The invariant di-boson mass spectrum is shown in Fig. [3 where 
the Higgs resonance is clearly visible for all Higgs masses. 

In order to determine the resonance parameters (m, T, N) from this distribution, the theoret- 
ical Breit-Wigner shape of the resonance has to be convoluted with a properly tuned detector 
resolution function. The detector resolution is estimated by using MC with zero Higgs width, 
Th = and the same selection as described before. 
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Figure 5: Di-boson mass resolution from MC with Fh = OGeV. As can be seen, resolutions 
differ for H WW and H TIL. Both distributions are asymmetric due to the constraints of 
the kinematic fit. 



It turns out that this detector resolution is different for H — > WW and H 
which enter the di-boson mass spectrum with relative fractions 



ZZ decays, both of 



RFh- 



A'h^ww 



►WW 



A'h^ww + Ah_»zz 



and 



RFu- 



»zz 



A^H-^ww + Ah^zz 



respectively, with A^ being the number of events for each channel after selection. The different 
resolutions arise from the fact that for H — > ZZ the correct di-boson mass partially is an iiqq 



6 



invariant mass, while for H WW it is always qqqq. The different expected mass spectra and 
the parameterizations^ used for convolution are shown in Fig. [5] for mn = 200 GeV. Resolutions 
for the other Higgs masses are similar. 

The fractions iJFn^ww ^-nd RF^i^zz are determined from the di-jet mass as obtained from the 
kinematic fit. An example for mn = 200 GeV is shown in Fig. El where two peaks from W- and 
Z-decays respectively are clearly visible. However, the tails are overlapping. 
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Figure 6: Jet-jet mass as obtained from the kinematic fit. Clearly visible are two peaks at mw 
and mz respectively. 



The di-jet mass spectrum is divided by a (in principle arbitrary) mass cut, choosen to be mjj = 
86 GeV. The number of events from all event classes (H — > WW, H — > ZZ and background), Ni 
is broken up in events below and events above this cut. The corresponding probabilities 
are referred to as = N^/Ni. Expectations for and N^kgy^ are derived from MC, while 
the total number are counted in data. 

From this, the relative fractions can be calculated according to: 
= ^bkgrd + ^H^WW ^l|_WW + ^H^ZZ -Pr-^zZ 

= Kkgrd + (^H^WW + iVH^Zz) {RFn^wW Pr^wW + -R^H^ZZ ^jt^zz) 

= Kkgrd + {Nu^WW + Nn^zz) ((1 " RFu^zz) P^^WW + ^^H-.ZZ P^^zz) ■ 



With iVn^ww + A^'h^zz = Ntot - N^kgrd resulting in 

Kot = Kkgrd + (Ntot - Ntkgrd) ((1 " iJ^H-^Zz) P^^^^ + RFu-.ZZ P^^zz) 
j^p^ _ ^ioi ~ ^bkgrd -Ph_^WW 



and iJFn^ww = 1 — RFh 



{Ntot Nbkgrd) (-fll-^ZZ -^H-^Ww) -^H-^ZZ -^H-^WW 

^tot ~ ^bkgrd Pu^ZZ 



+ZZ 



{Ntot Nbkgrd) {Pr^ww Pr^zz) Pr^ww Pr^zz 



^Detector effects are parameterized separately for H — > WW and H ^ ZZ by multi-gaussian functions. The 
choice is arbitrary and motivated by the good agreement between MC and parameterization. 
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Due to the larger background contributions below the cut, the determination of RF-n^zz and 
calculation of RFn^^^y^i = 1 — RFu^zz is more accurate than the opposite way. Nevertheless 
the analogue determination of -R-Fh^ww by considering events below the threshold can be used 
as a cross-check. 

The determination of the relative fractions is model independent. In addition, assuming there 
are only Higgs decays to W- or Z-pairs, the corresponding branching rations BRh^ww and 
BRh-+zz can be calculated. For the calculation, the selection efficiencies en^ww and en^zz are 
taken from MC. 

The number of events selected is a function of cross section, integrated luminosity, selection 
efficiency, and branching ratios. Consequently, the relative fraction iiFn^zz can be expressed 
as 



RFu^ 



►zz 



o"HZ en-tzz BRh^zz BRzzz_,2£4q 



BRzWW^2^4q) 



o"HZ (eh^zz BRh^ZZ BRzzZ^2£4q + ch^ww BRh^ww 
_ en^zz BRh->zz BRzzz^2^4q 

en^zz BRh^zz BRzzz^2£4q + ^h^ww BRh^ww BRzww^2f4q 

With BRh^ww + BRh^zz = RFu-,ww + RFn^zz = 1 and BRzzz^2^4q = 3BRzww^2^4q this 
transforms to 



BRh- 



ch^ww RF}i-^zz 



>zz 



eH-»ww RFn-^zz + 3 en^zz RFr^ww 



The factor 3 arises from the threefold ambiguity in ZZZ 2£4q. 

Errors for relative fractions and branching ratios are calculated using error propagation. All 
parameters derived from MC are assumed to be known exactly, so that AN^ 
only errors taken into account. Values obtained are listed in Tab. 01 
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mu [GeV] 


ARFh^ww 


RFh~.zz 


ABRh^ww 
BRh->ww 


ABRH-.ZZ 
BRh^zz 


200 


7.3% 


5.7% 


3.5% 


9.9% 


240 


8.9% 


7.0% 


5.0% 


10.8 % 


280 


11.9% 


9.3% 


7.7% 


16.2% 


320 


15.2 % 


12.0% 


8.6% 


17.3 % 



Table 3: Results for relative fractions and branching ratios 



5 Reconstruction of Higgs Resonance Parameters 

Finally, the spectrum of the reconstructed di-boson masses is built. Each event yields three 
entries for the different combinations possible. This spectrum can be described by following 
parts: 

1. The correct combination of di-bosons form a clear peak. This peak can be described as 
a theoretical Breit-Wigner function convoluted with the detector resolution. The Breit- 
Wigner parameters /i, T and are free parameters of the fit, while the parameters of the 
detector resolution are fixed from MC. 

2. The wrong combinations of di-boson masses form a flat combinatoric background. The 
shape is parameterized by a step function whose parameters are fixed from the same MC 



8 



sample as the detector resolution^ while the number of entries is determined in the fit by 
the free parameter of the Higgs peak Breit-Wigner function. 



3. In addition, there is a flat distribution from background events left after the selection. 
This physical background is parameterized by a similar step function as the combinatorical 
background. In this case besides the parameters for the shape also the number of entries 
is fixed by MC expectation^. 

Fig. [7| shows the distributions of reconstructed di-boson masses with the parts described and the 
fitted function. The fit parameters /i = mn, T = Fh and N = a x BR x e are determined in 100 
independent MC experiments, each corresponding to f£ = 500 fb^^ of integrated luminosity. 
Both mean values and spread of the results lay within statistical expectations, mean values for 
the measurement precision are listed in Tab.0J 
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Figure 7: Di-boson mass spectra for the four Higgs masses under study. Everywhere, the higgs 
resonances (yellow) are clearly visible over physical (blue) and combinatorical (green) background 
for jC = 500 fb~^ of integrated luminosity at ^/s = 500 GeV. 



6 Summary and Conclusion 

We present a method for measuring Higgs mass, width and event rate in a model independent fit 
from the reconstructed Higgs lineshape at TESLA. The method is restricted to Higgs bosons with 

^Here it is assumed that the shape of the combinatorical background does not depend on the Higgs width. 
^In later experiments, the background parameters can be determined off-peak in data. 
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7.7% 


16.2% 


320 GeV 


6.3% 


0.36% 


26.4 % 


8.6% 


17.3 % 



Table 4: Resolutions for event rate N , Higgs mass ttih and Higgs width Fh as fitted to the di- 
boson mass spectrum. Also listed are results for the branching ratios BRh^ww o,nd BRh^zz is 
described in Sec. ^ All numbers are mean values obtained with 100 independent signal samples 
corresponding to JC = 500 fb^^ at y/s = 500 GeV each. 

widths as large as few GeV. Otherwise, precision is spoiled by detector mass resolution which 
is in the same order. Assuming the Higgs decays to W- and Z-boson pairs only, determination 
of the corresponding branching ratios BRh^ww and BRh^zz is possible as well. 

The selection is based on reconstructing e'''e~ — > HZ events with successive H WW/ZZ 
decays. Selecting final states with two charged leptons and four jets gives a good handle on 
background suppression and high event rates. Since identification of r- leptons is not modelled 
in the fast detector simulation used, only electrons and muons are considered. 

Four jets are paired to two boson candidates by a 4C kinematic fit, two identified leptons form a 
third boson candidate. The Higgs resonance is reconstructed as di-boson mass. Higgs event rate, 
mass and width are extracted from the resonance lineshape in a model independent fit. The 
results of 100 independent MC experiments each corresponding to /£ = 500 fb^^ of integrated 
luminosity lay within statistical expectation. Mean values of the precision achieved are listed in 
Tab.H 

Main motivation for this study was to explore a direct method for measuring the total Higgs 
width Fh. Before, this has only be studied for LHC experiments jlSi- Results obtained in this 
study are comparable in precision for mass and width of the resonance. No LHC numbers are 
available for event rates, and branching ratio determination is impossible from the lineshape 
reconstruction. 

Fig. El (left) compares the results obtained in this study on Higgs width and mass with previous 
TESLA studies on lower Higgs masses. For the Higgs width, the indirect determination via 
the Higgs coupling to W-bosons, measured in the cross section of WW- fusion jjj, is the most 
precise method studied so far. However, up to now only H ^ bb decays have been investigated, 
so precision breaks down as does the branching ratio BR^^i^^ at mn ^150 GeV. On the other 
hand, precision for direct width determinations in the Higgs lineshape are limited by the narrow 
Higgs width below mn^: 200 GeV. The gap could be closed by indirect determinations with the 
analysis of H ^ WW decays. The hope is, that combination of direct and indirect measurements 
significantly improves the precision on the Higgs width for high Higgs masses. 

Precisions of the direct method alone may be optimized as well. There are many ways to 
enhance signal event rates by considering more common final states. For instance, one could 
include Z — > r+r" or study HZ qqqqqq and HZ — > qqqq^i/ final states. But, the effects of 
larger background contributions and final state neutrinos need dedicated studies. 

It is pointed out that this study is optimized for width determination. Results on mass, event 
rate and branching ratios are welcome spin-offs, but there might be other processes and methods 
still to be studied which are more appropriate. As an example, the precision on Higgs mass 
measurements of previous TESLA studies Jl] is compared to results of this study in Fig. |H1 
(right). As can be seen, precision is about three times better in the dedicated study. The main 
reason, is selection of more common signal final states and thus higher statistics. 
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Figure 8: Precisions of Higgs width (left) and mass (right) measurements at TESLA. Results 
from former studies are shown in red, those of this simulation in blue. 
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